Transmission and Scattering of a Lorentz Gas on a Slab 
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We perform numerical scattering experiments on a Lorentz 
array of disks centered on a triangular lattice with L columns 
and study its transmission and reflection properties. In the 
finite horizon case, the motion of the particles may be mod- 
eled as simple one dimensional random walks with absorbing 
walls for which the scaling of the transmission and reflection 
coefficients are known, and agree with those found numeri- 
cally. In the infinite horizon case the analogy with a simple 
diffusive process is no longer valid. In this case we compare 
our results to those expected for a one dimensional Levy walk, 
again with absorbing walls, for which logarithmic corrections 
to the scaling relations appear. These corrections are con- 
sistent with the numerical results. The scaling with L and 
the symmetry properties of the forward ot{4>) and backward 
cm((f>) differential cross sections are also studied, and some of 
their salient features are discussed. 

PACS number(s): 05.40 +j,05.45 +b,05.60 +w 



I. INTRODUCTION 

The Lorentz gas is an ensemble of noninteracting point 
particles which move freely with clastic reflections from 
fixed scatterers M. It is a basic model for linearized 
kinetic equations |2| and its ergodic properties are well 
known ||. 

In this paper we present the results of numerical exper- 
iments in which a large number of particles are incident 
on an array of disks centered on a triangular lattice. The 
particles are launched initially either in the +x direction 
or isotropically towards the array, and are reflected elas- 
tically from the scatterers. The array of disks is finite 
in the x direction and infinite in the y direction so we 
speak of a "slab" of scatterers. Our slab is characterized 
by two parameters, the width of the slab, that is, the 
number L of columns of scatterers, and the separation w 
between them (the disk radius is unity). The quantities 
that are measured are the transmission T and reflection 
R coefficients, the mean survival time r of particles in the 
slab and the transmitted ot and reflected <tr differential 
cross sections. These quantities are analyzed as functions 
of the parameters characterizing the slab. 

If the separation between scatterers w in a triangu- 



lar lattice is small, w < w c = 0.3094 . . ., the length of 
free motion of the point particles is bounded, that is, the 
particles "see" a finite horizon. On the other hand, when 
w c < w the length of free motion may be unbounded, the 
particles "see" an infinite horizon. In this work we study 
the scattering properties in both situations, and analyze 
our results in terms of the characteristic motion of the 
particles in each case. In the finite horizon case, the mo- 
tion of the particles is known to be diffusive Q] and the 
diffusion coefficient can be estimated with ergodic argu- 
ments Q . Velocity correlation functions have been shown 
to decay exponentially [Q , as confirmed by numerical ex- 
periments |5J. In contrast, for the infinite horizon case 
the diffusion coefficient diverges logarithmically [j^U and 
correlation functions have a power law decay || . We find 
that there are also fundamental differences in the scat- 
tering properties for each case. 

Lorentz gases in finite size geometries have also been 
introduced elsewhere with the aim of studying escape 
rates and their relation to transport coefficients and frac- 
tal repellers An account of these results, together 
with a formulation of the problem in terms of flux bound- 
ary conditions can be found in Ref . . 

In Section |l| we introduce all the definitions and dis- 
cuss in detail the numerical experiments for the case of 
a finite horizon. Section III draws an analogy between 



the Lorentz scattering experiments in the finite horizon 
case and the behavior of one-dimensional diffusive motion 
with absorbing boundaries. This serves as a basis for the 
explanation of the observed scaling laws for transmission 
coefficients and survival times. In Section |^ we present 
the results for the infinite horizon case when the particles 
are launched initially in the x direction. These show log- 
arithmic corrections to the scaling laws, which are consis- 
tent with considering the motion of the particles within 
the slab as a Levy walk. In Section |v| we discuss angu- 
lar dependences and symmetries of the transmission and 
reflection differential cross sections. In Section VI we dis- 



cuss how some of the scattering properties are affected by 
sending the particles isotropically, i.e. with an incidence 
angle uniformly distributed between —it/ 2 and tt/2. A 
modified random walk model displaying these same dif- 
ferences is also briefly discussed. Finally, Section |VIl| is 
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devoted to discussion of the results and conclusions. 



II. SCATTERING WITH A FINITE HORIZON 

The scatterers we consider are disks of unitary radius 
centered on a triangular lattice as shown in Fig. |l|. The 
distance between neighboring centers is 2 + w and the 
centers lie along L lines parallel to the y axis. The slab 
is finite in the x direction and infinite in the y direction. 
A large number N of particles are incident from the left 
parallel to the x axis with unit speed. The particles move 
freely except for elastic collisions at the boundary of the 
disks. In the experiments, the dynamics is solved by 
considering the motion in the elementary Wigner-Seitz 
hexagonal cell where opposite sides are identified. Each 
incident particle has a different impact parameter b, de- 
fined here as the distance between the initial position 
and the horizontal line passing through the center of the 
scatterer in the Wigner-Seitz cell. Due to the symme- 
try of the slab, it is sufficient to consider b between 
and 1 + w/2. The trajectory in the slab is obtained by 
unfolding the orbit in the Wigner-Seitz cell. The cases 
where particles are incident with an angle different from 
zero and the effect of isotropic incidence will be briefly 



discussed in Section VI 
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FIG. 1. Slabs of scatterers in a triangular array, (a) Finite 
horizon, w — 0.3, with L = 5. (b) Infinite horizon, w = 1., 
with L = 5. 

In the finite horizon case, < w < w c = (4/V3- 2) = 
0.3094 . . ., the particles that enter the scatterer cannot 
travel long distances without suffering collisions. In the 
infinite horizon case, w c < w, the particles may travel 
arbitrarily far between collisions, due to the opening of 
infinite corridors. If w c < w < 2 the infinite corridors 
lie at angles of 7r/6, 7r/3, and tt/2. We will be mainly 
considering particles incident parallel to the x axis that 
cannot cross the slab without collisions when w < 2. 

Since the slab is infinite in the y direction and the col- 
lisions are elastic, every particle that enters the slab must 
eventually exit it, except for a set of measure zero which 
goes asymptotically to periodic orbits inside the slab (we 



disregard all zero measure sets herafter). Thus, in prac- 
tice, a particle that enters the slab collides with some of 
the obstacles and will be ultimately transmitted or re- 
flected, leaving the slab with an angle <f> measured with 
respect to the +x direction. Particles are transmitted if 
they exit the slab from the right (\<fi\ < tt/2), and are 
reflected if they exit from the left (\<p\ > tt/2). 

A first characterization of the system is through the 
computation of the transmision T, and reflection R, coef- 
ficients. The former is defined as the fraction of particles 
that pass through the slab and exit on the right, and the 
latter as the particle fraction that exits the slab on the 
left (obviously T + R = 1). A finer quantity is the dif- 
ferential scattering cross section a defined by saying that 
a(4>)d(f) is the fraction of particles scattered between (f> 
and 4> + d<j). We can separate this quantity in the trans- 
mitted <jt and reflected ctr differential scattering cross 
sections by considering \<j>\ < 7r/2 and \<j)\ > tt/2 respec- 
tively. 

In Fig. ^ we show the dependence of the transmission 
coefficient T on L for the finite horizon case. These re- 
sults are consistent with the scaling law T ~ L" 13 where 
P « 1 independently of w. This behavior will be justified 
in the next section drawing from an analogy with random 
walks. Since R = 1 — T and T scales to zero with L, R 
does not depend on L for L > 1. Another quantity of 
interest is the mean survival time inside the slab r. This 
average time is evaluated over all the N incident particles 
no matter if they are transmitted or reflected. We find 
that r ~ L 7 with 7 w 1 as we show in Fig. a. 



III. DIFFUSIVE BEHAVIOR IN FINITE SIZE 
SYSTEMS 

In the finite horizon case, rigorous results, convincing 
evidence and plausible arguments have been set forth in- 
dicating that the motion of the particles in the Lorentz 
system can be accurately modeled as a simple random 
walk . In its simplest version the particles can be 
viewed as hopping between adjacent "cages" in an essen- 
tially uncorrelated fashion, and staying in each cage for 
a well defined average time. For the case under study 
in this paper it is not even necessary to consider the 
random walk process as occuring on a two dimensional 
lattice since the quantities we are interested in can be 
calculated from the projection of the walk onto the finite 
direction x. 

While the discrete one dimensional random walk on 
a finite lattice can be described completely |lj,[l3| , such 
a detailed comparison between the two systems cannot 
hold. As the random walk is an analogy to the Lorentz 
system, the best we can realistically expect to deter- 
mine from it is the scaling behavior of the quantities 
under study. With this in mind, we choose to evaluate 
the transmission and reflection coefficients for the simple 



2 



random walk problem via the diffusion equation, with a 
numerically determined phenomenological diffusion con- 
stant D. This equation describes the evolution of the 
coarse grained particle density in the system and is ex- 
pected to be valid when the system size is much greater 
than the root mean square (rms) step length. Once again, 
since the slab is translationally invariant along the y axis, 
the coarse grained particle density obeys a diffusion equa- 
tion along x. 
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FIG. 2. Logarithmic plot of the transmission coeffi- 
cient T as a function of the size of the system L, for dif- 
ferent values of w. The number N of incident particles 
were, for w = 0.1,4 x 10 7 ; for w = 0.15,2 x 10 7 , and for 
w — 0.2, 0.25, 0.3, 10 7 . The full lines are the least square 
fits, all of them compatible with the theoretical prediction 
T ~ L' 1 to within 1%. 
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FIG. 3. Logarithmic plot of the survival time r vs. L for 
the same values of w and N as in Fig. [| The full lines are the 
least square fits, all of them compatible with the theoretical 
prediction r ~ L to within 1.6%. 

To estimate the reflection and transmission coefficients 
within this approximation, we require the solution of the 
diffusion equation with absorbing boundary conditions at 
x = and x = X and a constant unit input flux at site a. 



In the steady state, the magnitude of the fluxes at and 
X will give the splitting probabilities; i.e. the probability 
that a particle injected at position a will be absorbed at 
the origin or at X. In our scattering system the particles 
are inciding on the left, which can be thought of as having 
the injection point a near the origin. Then the calculation 
outlined above yields the transmission coefficient 

T = a/X. (1) 

To estimate the mean survival time of the particles 
within the slab, we recall that within the diffusion ap- 
proximation, the survival time for a random walker start- 
ing at position a satisfies the equation 

1 (2) 

with the conditions r(0) = t(X) — 0, where D is again 
the diffusion constant of the process. Thus 

r(a) = ^a(X-a) , (3) 

and if the injection point is taken to be close to the origin 
(o<I),we obtain 

r - aX/2D. (4) 

These results can also be obtained on a more general 
footing through Wald's identity |ll|, and are expected 
to hold as long as the rms step size of the random walk 
is small compared to the system size X , and the number 
of steps given by the random walker scales linearly with 
time (i.e. there are no long tail waiting time distribu- 
tions). 

We can identify the quantities appearing in Eqs. (|l|) 
and (^) corresponding to the Lorentz scattering ex- 
periment. The diffusion coefficient D is computed in 
Ref. fa] in a random walk approximation and numeri- 
cally through the Green-Kubo relation. The length X is 
related to w and L by X — L(2 + w)V3/2 and thus the 
penetration length a can be determined by the slope of 
the dependence of r on L as in Fig. |^. The results are 
consistent with the distance between traps (2 + w)/Va, 
as defined in Ref. 0] , for small values of w (where also the 
diffusion constant predicted in the random walk approx- 
imation agrees with the numerically determined one) . 

IV. INFINITE HORIZON 

When the horizon becomes infinite the analogy to the 
simple diffusive process breaks down. This occurs as a 
consequence of the opening of infinite corridors between 
scatterers in which the particle is capable of travelling 
very large distances between collisions. The distribution 
of the length of these sojourns, p(r), has been shown to 
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oo by both numerical results 



decay as r~ 3 when r 

and theoretical arguments J? 14 |Jl^]. Thus the rms step 
length diverges and the diffusive approximation described 
in the last section breaks down. 

If we insist on making a random walk description of 
the system, we are now led to consider a random walk 
with a distribution of step lengths without second mo- 
ment (generically called "Levy nights" [ ^3p^ ]). Here 
one sometimes makes the distinction between a discrete 
time process, in which each step is selected from a power- 
law distribution of lengths, but always takes a fixed time 
and the so-called Levy walk in which each step takes a 
time proportional to its length. This latter, while more 
relevant to the case we are discussing, is not significantly 
different from the Levy flight if the first moment of the 
step size distribution exists, as it certainly does in our 
case. We shall therefore ignore the considerable com- 
plications this causes and often identify time with the 
number of jumps or collisions n. 

In contrast to the diffusive case, the derivation of the 
transmission coefficient in the case of Levy flights appears 
not to have been treated in the literature. We present 
an argument which leads to a scaling prediction which 
seems reasonable for general Levy flights, and specialize 
it to the case we are concerned with. 

Denote the step distribution of the Levy flight as 
P(r) ~ l/|r| 1+Q with < a < 2. A random walker 
with this step distribution travels a typical distance Jl6| 

f t 2 / a for 1< a < 2 
£ 2 ~ I tint for a = 2 
W for a > 2 

in time t. 

If we assume that the particle is initially at a point a 
sufficiently near the origin (a <C A), then we can invoke 
Sparre- Andersen's theorem |l2|] . This theorem states 
that the probability that the walker steps for the first 
time to the left of the point a at the n-th step is a uni- 
versal function 7r(n), which is completely independent of 
any of the properties of P(r), as long as P{r) is symmet- 
ric and continuous. This distribution is found to decay 
as 7r(n) ~ -nT^I 2 for any kind of unbiased walk whatso- 
ever, in particular for the Levy flights we are concerned 
with. This allows us to estimate straightforwardly the be- 
haviour of Levy flights both as regards their transmission 
and the mean survival time in the interval considered. 

The scaling behavior of the transmission coefficient of 
a Levy walk across a finite system of length X can be esti- 
mated as follows: Denote by r a (X) the time required for 
the walker to travel a distance of order X with apprecia- 
ble probability. The transmission coefficient is then the 
probability that the walker never steps left of the origin 
during r a (X). From the above scaling relation, T a (X) is 
expected to scale as X a for 1 < a < 2, as A 2 / In A for 
a = 2 and as X 2 for a > 2, which corresponds to nor- 
mal diffusion. Then, in terms of t q (A) the transmission 



coeficient is given by 
T(X) ~ V 7r(n) ~ / t^dt ~ l/^JX). 



n>T a (X) 



(5) 



For a > 2 we therefore get the behavior of the or- 
dinary diffusive case treated in Section III. In partic- 
ular, for our infinite horizon Lorentz slab, we expect 
T(X) ~ y/hTZ/L. 

As for the survival time, it is estimated in a similar 
way: The probability of leaving the interval [0, X] at time 
n is given by ir(n) as long as n is not so large that leaving 
at the right hand-side becomes appreciably likely. This 
occurs at times of the order of t q (A). This reasoning 
gives for the mean first exit time r 



r a (X) 
n=l 



Vr a (X). 



(6) 



The prediction for the probability distribution 7r(n) of 
leaving the slab to the left after n collisions based on the 
Sparre- Andersen theorem can be tested for the Lorentz 
gas; it is valid both for the finite horizon case and for 
the infinite horizon. In Fig. ^| we show this distribution 
as obtained from the simulation for w = 1.5 (infinite 
horizon) together with the Sparre- Andersen scaling law. 
The agreement is good even for small values of n. 




FIG. 4. The probability distribution 7r(n) of leaving the 
system to the left after n collisions, for w = 1.5, and from 
left to right L = 100, 500, 1000, 1500 and TV = 10 7 . The 
Sparre- Andersen scaling result is shown by the dashed lines. 

The scaling laws for the transmission coefficient and for 
the mean survival time can be also tested for the Lorentz 
gas in the infinite horizon case. In Fig. | we show (TL) 2 
for a set of w values as a function of ln(L). According to 
our theory {TL) 2 ~ ln(L) for w c < w. In Fig.||we show 
(t/L) 2 as a function of ln(L) in order put in evidence the 
predicted logarithmic correction. 
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Runs were also carried out at other fixed incidence an- 
gles, and the same features as described above were ob- 
tained. We have thus shown that, given a fixed angle 
incidence, the opening of the horizon appears to produce 
logarithmic corrections to the scaling laws present for fi- 
nite horizon. This feature is shared by the behavior of 
other quantities which also present logarithmic correc- 
tions, such as the diffusion coefficient. 
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FIG. 5. The logarithmic correction to T is shown by plot- 
ting \{LT{L)f - (100T(100)) 2 ]/(100T(100)) 2 vs. ln(L) for 
w = 0.1 (with N = 4 x 10 7 ), and w = 0.6, 1.0, 1.5, 2.0 (with 
N — 10 7 ). Straight lines with nonzero slope, indicating the 
logarithmic correction, are expected for w above the infinite 
horizon. The dotted lines are meant as a guide to the eye. 
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FIG. 6. The logarithmic correction to r is shown by plot- 
ting \{L/t{L)) 2 ) - (100/t(100)) 2 ]/(100/t(100)) 2 vs ln(L) for 
the same values of w and N as in Fig. H. Again straight lines 
with nonzero slope are expected for w above the infinite hori- 
zon, and we also show the corresponding quantity por a Levy 
walk with a = 2.0. The dotted lines are meant as a guide to 
the eye. 



V. ANGULAR DEPENDENCE OF FORWARD 
AND BACKWARD SCATTERING 

We have already observed that for fixed incidence an- 
gle the transmission coefficient scales as 1/L in the fi- 
nite horizon case, whereas it scales as y/hiL/L in the 
case of infinite horizon. It is also instructive to look in 
somewhat greater detail at the angular distribution of the 
transmitted and reflected particles (this is what is done 
also in chaotic scattering experiments with fewer scat- 
terers 0). We define a R {4>) as the density of particles 
reflected at angle ef> and ut{<P) similarly for the transmit- 
ted particles. Due to the symmetry 0^—0 the range 
of <f) will be from to ir/2 for transmitted particles and 
from 7r/2 to ir for those that are reflected. As these dis- 
tributions are also dependent on L, we further propose 
that in the finite horizon case they can be expanded as 



(Tr(4>) = a R ((j)) 



M</>) 



and 



a T {4>) = 



b T {<l>) 
L 



(7) 



(8) 



The leading term in the transmission cross section is zero 
since there is no transmission for infinite L. The follow- 
ing relation has then been found to hold in all cases for 
sufficiently large L (see Fig. 0): 



M0) = -frrO - 4>) 



(9) 



This symmetry relation can be understood as follows: 
Imagine that the slab is subjected to a continuous flow 
of particles with identical distributions incident from 
both sides of the sample. If we now look at the distri- 
bution of angles of the particles which go to the right 
in these circumstances, it is given by the superposition 
arty) + o- R (4>) = a R {<j>) + 4- M + • • , Now, by 
symmetry, the above set-up is equivalent to having par- 
ticles incident only from the right and a reflecting wall 
at the middle of the sample. However, from the fact that 
a particle travelling deep into the slab eventually loses 
memory of its original incidence direction, it follows that 
the trajectory of a particle reflected on the wall is indis- 
tinguishable with the trajectory of a particle travelling in 
a semi-infinite system and eventually returning across the 
position of the wall (which occurs with probability one). 
Thus, given that the correlation with the initial incidence 
is small enough, the system with the reflecting wall will 
not show any great difference from the semi-infinite sys- 
tem (L = oo) as far as the angular distribution of its 
particles is concerned. This then allows to derive the 
identity 



a R ((j)) = a R {(j>) 



b R (<j>) + b T (ir- <j)) 



(10) 
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from which the result follows. Numerically, this is quite 
well borne out both in the case of finite (Fig. [?]) and 
infinite horizon (Fig. . 
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FIG. 7. The quantities b T ((p/Tv) (•) and 6r((tt - 0)/tt) 
(A) vs. 0/tt for w = 0.3 and TV = 10 7 . The values for 
br are extracted from data with L — 100 and those for 6_r 
are determined by substraction of data with L = 200 and 
L = 100. 




FIG. 8. The quantities -b T (cf>/ir) (•) and b R ((ir - 4>)/n) 
(A) vs. 4>/tv for w — 1.5 and TV = 10 7 . The values for 
br are extracted from data with L — 100 and those for ba 
are determined by substraction of data with L — 200 and 
L = 100. 

In the case of infinite horizon, another striking feature 
of the angular distribution function is found: there is a 
clear peak in bx around the value 4> c = tt/6 as well as a 
corresponding dip in an around the value <f> c = 5ir/6 
(see Fig.^) which are the angles the infinite corridors 
make with the edge of the sample. Qualitatively, the ap- 
pearence of these singularities can be argued as follows. 
As is well-known, the periodic system corresponding to 
the one we are studying is ergodic, so that all allowed po- 
sitions and all directions are eventually equally probable. 
In the finite system we are studying, this is no longer the 



case. In particular, near the edges, the relative weight 
of the directions leading to escape will be different from 
the other ones. Nevertheless, we may at first start with 
the approximation that, at least as long as the particles 
enter reasonably deep into the system, the hypothesis 
of cquidistribution of velocity directions holds to a good 
approximation. Then we may estimate the transmission 
coefficient as a function of angle by means of the fraction 
of surface area from which trajectories escape at that 
angle (we may assume that the transverse direction of 
the gas is made finite by some device such as periodic 
boundary conditions). This yields a singularity in the 
angular dependence of the transmission coefficient near 
the critical angles when L — > oo. 




FIG. 9. The quantities br (full line) and an (dashed line) 
vs. 4>/tv for w — 1.5, TV = 10 6 . The features expected at 
4> — 7r/6 and <f) = 5n/6 are clear. There is also a peak at 
4> ~ arctan\/2, for which we have no explanation. 



VI. THE EFFECT OF ISOTROPIC INCIDENCE 

The following apparent difficulty motivated us to also 
study the case in which the particles are isotropically inci- 
dent upon the slab: In the preceding Section, we reported 
numerical and theoretical evidence for the existence of a 
logarithmic correction to the mean survival time. Yet a 
quite general argument appears to show such a correction 
is impossible. Consider the phase space on a constant en- 
ergy shell inside the Lorentz array. The volume of this 
phase space, which clearly scales as L, can be expressed 
as the integral of the time of residence inside the sys- 
tem over all points of entry (this relation is known as 
the Katz formula jL8|). Since the volume of trajectories 
that remain forever confined to the system is zero, and 
the volumes of those entering from the left and from the 
right are equal, it follows that the mean survival time 
scales as L, in contradiction to what was numerically ob- 
served in the last Section, as well as to the predictions of 
the Levy flight model. 

If we consider the above argument carefully, however, 
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we see that it only applies if all angles of incidence are 
taken as equally probable. To understand why this might 
make a difference, one should note the following: If all an- 
gles are equally probable, then the distribution of initial 
step lengths will have a singularity due to the probability 
of launching the particle with an angle very close to crit- 
ical and an appropiatc impact parameter. This leads, as 
is readily seen, to the probability for a large initial step 
length of x going as x~ 2 . This behaviour is in marked 
contrast to the step probability distribution within the 
interior of the system, for which the impact parameter 
and the angle are interrelated. For these, as was pointed 
out before, the probability of a large step x goes as a; -3 . 

In order to clarify these issues, we performed simula- 
tions on the system with isotropic incidence. Indeed, as 
expected from the above argument, the mean survival 
time was found to scale as L, without any logarithmic 
corrections (see Fig. ^). On the other hand, the transmis- 
sion coefficient was found to retain its peculiar behaviour 
(see Fig. |^). The above argument is therefore sound, it 
does not contradict the numerical work reported in the 
previous Section. On the contrary, since the two mod- 
els show such clear differences with regard to their mean 
survival time, this seems to indicate that the distribution 
of initial step lengths may well have been the cause. 

To test this final hypothesis, we also simulated a Levy 
walk with a = 2 in which the first step has a distribution 
with an x~ 2 tail, corresponding to a = 1. The results 
are again very clear: the mean survival time now grows 
as L without any logarithmic corrections. On the other 
hand, the transmission coefficient still has the previous 
anomalous behaviour, though it takes a somewhat longer 
time to reach it. 

VII. CONCLUDING REMARKS 

In this work we have examined the scattering proper- 
ties of a Lorentz gas incident on an array of scatterers 
centered on a triangular lattice with a finite number L 
of columns. Much of our attention has focused on the 
scaling with L, for large values of L, of transport and op- 
tical properties, namely, transmission and reflection coe- 
ficients, mean survival time and differential cross section. 
It should be emphasized that though most of the numer- 
ical results reported in this paper were obtained for inci- 
dence along the x direction, exactly the same phenomena 
was observed in runs carried out at other fixed incident 
angles. On the other hand, some significant differences 
were found when the particles were launched isotropically 
in the infinite horizon case. 

For our understanding of the observed numerical 
trends we have considerably profited from the link, in 
some cases formal, in others qualitative, to random walk 
processes. Two regimes are considered: the finite and the 



infinite horizon cases. For finite horizon the transmis- 
sion coefficient scales with 1/L, the survival time with 
L and the differential cross-section has no singularities 
and presents certain symmetry properties. All this is in 
agreement with the behavior of normal diffusion and or- 
dinary random walks with absorbing boundaries. 

The infinite horizon case, on the other hand, exhibits 
logarithmic corrections to the aforementioned scalings 
depending on the incidence angle distribution (in this 
case the relation to Levy walks is illuminating). Also, 
singularities appear in the differential cross section at 
angles corresponding to the corridors inside the slab, for 
which we only have a partial understanding. A peculiar 
effect related to the difference observed between particles 
launched in one fixed direction and particles launched 
isotropically can also be explained in terms of a simple 
Levy flight model with a different distribution for the 
initial step length. This allows an explanation for the a 
discrepancy in the behavior of the mean survival time, 
which had been anticipated on quite general grounds. 

It is also of interest to understand the features of the 
free motion length distribution. For finite horizon two 
peaks in such distribution are present at the values w and 
y/3(2 + w) — 2, which correspond to the unstable periodic 
orbits perpendicular to the disks. For infinite horizon a 
set of peaks develops, whose number increases with w; 
the slope of the envelope of the probability distribution 
is —3. The origin of such peaks, which most probably are 
related to other periodic orbits, remains to be analyzed. 
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